A new computational method for the detection of horizontal gene transfer events

In recent years, the increase in the amounts of available genomic data has made it easier to appreciate the extent by which organisms increase their genetic diversity through horizontally transferred genetic material. Such transfers have the potential to give rise to extremely dynamic genomes where a significant proportion of their coding DNA has been contributed by external sources. Because of the impact of these horizontal transfers on the ecological and pathogenic character of the recipient organisms, methods are continuously sought that are able to computationally determine which of the genes of a given genome are products of transfer events. In this paper, we introduce and discuss a novel computational method for identifying horizontal transfers that relies on a gene's nucleotide composition and obviates the need for knowledge of codon boundaries. In addition to being applicable to individual genes, the method can be easily extended to the case of clusters of horizontally transferred genes. With the help of an extensive and carefully designed set of experiments on 123 archaeal and bacterial genomes, we demonstrate that the new method exhibits significant improvement in sensitivity when compared to previously published approaches. In fact, it achieves an average relative improvement across genomes of between 11 and 41% compared to the Codon Adaptation Index method in distinguishing native from foreign genes. Our method's horizontal gene transfer predictions for 123 microbial genomes are available online at

A sensitive, support-vector-machine method for the detection of horizontal gene transfers in viral, archaeal and bacterial genomes

In earlier work, we introduced and discussed a generalized computational framework for identifying horizontal transfers. This framework relied on a gene's nucleotide composition, obviated the need for knowledge of codon boundaries and database searches, and was shown to perform very well across a wide range of archaeal and bacterial genomes when compared with previously published approaches, such as Codon Adaptation Index and C + G content. Nonetheless, two considerations remained outstanding: we wanted to further increase the sensitivity of detecting horizontal transfers and also to be able to apply the method to increasingly smaller genomes. In the discussion that follows, we present such a method, Wn-SVM, and show that it exhibits a very significant improvement in sensitivity compared with earlier approaches. Wn-SVM uses a one-class support-vector machine and can learn using rather small training sets. This property makes Wn-SVM particularly suitable for studying small-size genomes, similar to those of viruses, as well as the typically larger archaeal and bacterial genomes. We show experimentally that the new method results in a superior performance across a wide range of organisms and that it improves even upon our own earlier method by an average of 10% across all examined genomes. As a small-genome case study, we analyze the genome of the human cytomegalovirus and demonstrate that Wn-SVM correctly identifies regions that are known to be conserved and prototypical of all beta-herpesvirinae, regions that are known to have been acquired horizontally from the human host and, finally, regions that had not up to now been suspected to be horizontally transferred. Atypical region predictions for many eukaryotic viruses, including the α-, β- and γ-herpesvirinae, and 123 archaeal and bacterial genomes, have been made available online at

OMiR: Identification of associations between OMIM diseases and microRNAs

AbstractA large number of loci for genetic diseases have been mapped on the human genome and a group of hereditary diseases among them have thus far proven unsuccessful to clone. It is conceivable that such "unclonable" diseases are not linked to abnormalities of protein coding genes (PCGs), but of non-coding RNAs (ncRNAs). We developed a novel approach termed OMiR (OMIM and miRNAs), to test whether microRNAs (miRNAs) exhibit any associations with mapped genetic diseases not yet associated with a PCG. We found that "orphan" genetic disease loci were proximal to miRNA loci more frequently than to loci for which the responsible protein coding gene is known, thus suggesting that miRNAs might be the elusive culprits. Our findings indicate that inclusion of miRNAs among the candidate genes to be considered could assist geneticists in their hunt for disease genes, particularly in the case of rare diseases

The Transcription Factor Zfx Regulates Peripheral T Cell Self-Renewal and Proliferation

Peripheral T lymphocytes share many functional properties with hematopoietic stem cells (HSCs), including long-term maintenance, quiescence, and latent proliferative potential. In addition, peripheral T cells retain the capacity for further differentiation into a variety of subsets, much like HSCs. While the similarities between T cells and HSC have long been hypothesized, the potential common genetic regulation of HSCs and T cells has not been widely explored. We have studied the T cell-intrinsic role of Zfx, a transcription factor specifically required for HSC maintenance. We report that T cell-specific deletion of Zfx caused age-dependent depletion of naïve peripheral T cells. Zfx-deficient T cells also failed to undergo homeostatic proliferation in a lymphopenic environment, and showed impaired antigen-specific expansion and memory response. In addition, the invariant natural killer T cell compartment was severely reduced. RNA-Seq analysis revealed that the most dysregulated genes in Zfx-deficient T cells were similar to those observed in Zfx-deficient HSC and B cells. These studies identify Zfx as an important regulator of peripheral T cell maintenance and expansion and highlight the common molecular basis of HSC and lymphocyte homeostasis

The pandemic dilemma

Pairwise overlaps of TAD boundaries. The pairwise overlaps of TAD boundaries are shown for all samples of this study, after calling boundaries using hicratio (all reads, d = 0500). Before TAD calling, the Hi-C matrices were either unprocessed (filtered) or corrected using iterative correction (IC) (resolution = 40 kb). (PDF 3847 kb

Self-supervised learning in non-small cell lung cancer discovers novel morphological clusters linked to patient outcome and molecular phenotypes

Histopathological images provide the definitive source of cancer diagnosis, containing information used by pathologists to identify and subclassify malignant disease, and to guide therapeutic choices. These images contain vast amounts of information, much of which is currently unavailable to human interpretation. Supervised deep learning approaches have been powerful for classification tasks, but they are inherently limited by the cost and quality of annotations. Therefore, we developed Histomorphological Phenotype Learning, an unsupervised methodology, which requires no annotations and operates via the self-discovery of discriminatory image features in small image tiles. Tiles are grouped into morphologically similar clusters which appear to represent recurrent modes of tumor growth emerging under natural selection. These clusters have distinct features which can be identified using orthogonal methods. Applied to lung cancer tissues, we show that they align closely with patient outcomes, with histopathologically recognised tumor types and growth patterns, and with transcriptomic measures of immunophenotype

Notch signaling regulates metabolic heterogeneity in glioblastoma stem cells.

Glioblastoma (GBM) stem cells (GSCs) reside in both hypoxic and vascular microenvironments within tumors. The molecular mechanisms that allow GSCs to occupy such contrasting niches are not understood. We used patient-derived GBM cultures to identify GSC subtypes with differential activation of Notch signaling, which co-exist in tumors but occupy distinct niches and match their metabolism accordingly. Multipotent GSCs with Notch pathway activation reside in perivascular niches, and are unable to entrain anaerobic glycolysis during hypoxia. In contrast, most CD133-expressing GSCs do not depend on canonical Notch signaling, populate tumors regardless of local vascularity and selectively utilize anaerobic glycolysis to expand in hypoxia. Ectopic activation of Notch signaling in CD133-expressing GSCs is sufficient to suppress anaerobic glycolysis and resistance to hypoxia. These findings demonstrate a novel role for Notch signaling in regulating GSC metabolism and suggest intratumoral GSC heterogeneity ensures metabolic adaptations to support tumor growth in diverse tumor microenvironments

Transcriptional evidence for the "Reverse Warburg Effect" in human breast cancer tumor stroma and metastasis: Similarities with oxidative stress, inflammation, Alzheimer's disease, and "Neuron-Glia Metabolic Coupling"

Caveolin-1 (-/-) null stromal cells are a novel genetic model for cancer-associated fibroblasts and myofibroblasts. Here, we used an unbiased informatics analysis of transcriptional gene profiling to show that Cav-1 (-/-) bone-marrow derived stromal cells bear a striking resemblance to the activated tumor stroma of human breast cancers. More specifically, the transcriptional profiles of Cav-1 (-/-) stromal cells were most closely related to the primary tumor stroma of breast cancer patients that had undergone lymph-node (LN) metastasis. This is consistent with previous morphological data demonstrating that a loss of stromal Cav-1 protein (by immuno-histochemical staining in the fibroblast compartment) is significantly associated with increased LN-metastasis. We also provide evidence that the tumor stroma of human breast cancers shows a transcriptional shift towards oxidative stress, DNA damage/repair, inflammation, hypoxia, and aerobic glycolysis, consistent with the "Reverse Warburg Effect". Finally, the tumor stroma of "metastasis-prone" breast cancer patients was most closely related to the transcriptional profiles derived from the brains of patients with Alzheimer's disease. This suggests that certain fundamental biological processes are common to both an activated tumor stroma and neuro-degenerative stress. These processes may include oxidative stress, NO over-production (peroxynitrite formation), inflammation, hypoxia, and mitochondrial dysfunction, which are thought to occur in Alzheimer's disease pathology. Thus, a loss of Cav-1 expression in cancer-associated myofibroblasts may be a protein biomarker for oxidative stress, aerobic glycolysis, and inflammation, driving the "Reverse Warburg Effect" in the tumor micro-environment and cancer cell metastasis